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Abstract 

The solar and baseline neutrino oscillation data suggest bimaximal neutrino mixing 
among the first two generations, and trimaximal mixing between all three neutrino 
flavors. It has been conjectured that this indicates the existence of an underlying sym- 
metry for the leptonic fermion mass textures. The experimentally measured quantities 
however, are associated to the latter indirectly and in a rather complicated way through 
the mixing matrices of the charged leptons and neutrinos. Motivated by these facts, 
we derive exact analytical expressions which directly link the charged lepton and neu- 
trino mass and mixing parameters to measured quantities and obtain constraints on 
the parameter space. We discuss deviations from Tri-Bi mixing matrices and present 
minimal extensions of the Harrison Perkins and Scott matrices capable of interpreting 
all neutrino data. 



1 Introduction 



Since the experimental confirmation of neutrino oscillations, there have been assiduous 
efforts to measure the exact mixing angles and their tiny masses 0- The last few years 
we are in a position to know up to a high accuracy the three neutrino mixing angles and 
the mass squared differences. Several phenomenological explorations have led to the 
conclusion that the so called Tri-Bi (TB) maximal mixing ^ is in remarkable agreement 
with the solar and atmospheric neutrino data. Indeed, the experimental range of the 
three angles lie in the range 

sin^ 612 ~ 0.312l°:°i^, sin2 ^23 ~ 0.466^°;°^^^ sin^ ^13 « 0.126+°:°! (1) 

while the values of TB-prediction are 5,^,0 respectively. 

Another interesting aspect is the fact that in the basis where the charged lepton 
mass matrix is diagonal, TB-mixing is independent of the neutrino mass eigenvalues, 
the symmetric neutrino mass elements need only to satisfy three simple relations [T] 

One may attempt to attribute the regularity of the data in the leptonic sector to the 
existence of some particular symmetry of a suitable theoretical model. However, the 
above picture is definitely different from the corresponding one in the quark sector thus, 
in view of accumulating experimental evidence during the last decade, the reconciliation 
of neutrino data with simple U{1) family symmetry models |8] is rather unlikely. It 
has been suggested for example, that this specific structure might originate from a 
discrete non-Abelian symmetry [9l[in]. A different point of view is taken however in [7] 
where the authors claim that TB-mixing might be completely accidental as they found 
that significant violations from TB-mixing may occur within the present experimental 
bounds, with 1 — 3 mixing in particular leading to substantial deviations. Several other 
suggestions including the introduction of discrete and unified theories have appeared 
in the literature f3]- [2T]. 

The modifications on the TB mixing suggested by several of these proposals are 
based on perturbative considerations of the original TB-mass textures and the mixing 
matrix. However, in order to consistently study the effects on the mixing and the 
tiny mass differences in the neutrino sector, a rather accurate approach is required. 
Further, a major issue the present days is the exact measurement of the ^13 angle 
which in the original TB-model was assumed to be exactly zero, while data do allow 
for a small deviation. In addition, several measurable effects depend crucially on the 
value of the ^13 angle. For example, the survival probability of the reactor neutrinos 
involves both ^13 and Amg^ [22j, and current bounds allow a small value for ^13. If this 
angle is non-zero indeed, its value is sensitive to any modification of the mass matrix; 
in this case, approximate results and perturbative expansions may not be adequate 
to reliably determine the measured parameters in terms of mass textures eventually 
dictated by some symmetry. The aim of this letter is to fill this gap. Assuming only 

^For reviews see [T]- [5]. 
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general hermitean mass squared matrices for the charge leptons and neutrinos, we will 
derive analytical results for the mixing and mass-squared differences. To this end, using 
well known theorems from the spectral theory of matrices, we express a general 3x3 
hermitean fermion mass matrix as a second degree polynomial of a unitary matrix. 
This way, we are able to disentangle the mass eigenstates which appear only in the 
coefficients of this expansion. We then express the neutrino mixing angles as functions 
of variables that parametrize the unitary matrices which diagonalize the charged lepton 
and neutrino mass matrices H respectively. This procedure gives enough flexibility 
to determine the experimentally allowed range of the parameters, and, at the same 
time seek for mass textures eventually dictated by some underlying symmetry. As 
an immediate benefit of this approach we get a non-zero ^13 angle emerging even in 
the minimal TB scenario, provided that at least one non-zero phase in the neutrino 
texture is assumed. In addition, we will present a second example where compatible 
neutrino textures arise, otherwise not accessible by perturbative treatment around the 
TB solution. 



2 Formulation 



As explained in the introduction, the TB-maximal mixing is compatible with the ob- 
served neutrino data. However, the very specific form of the mass matrices postulated 
in this approach can only be embedded in particular classes of unified theories and even 
less string derived models. Early and present endeavors in this direction for example 
involve the rather promising A4 symmetry as far as the neutrino sector is concerned. 
Thus, ^4 can be generated by elements S,T satisfying S"^ = = (ST)^ = 1 and these 
can be viewed [15] as a subgroup of the modular group which plays a fundamental 
role in string theory. Nevertheless, a straightforward application to the quark sector 
is not very satisfactory since it predicts unacceptably small quark mixing, while only 
contrived variants can possibly reconcile data from both the quark and lepton sec- 
tors |16fl20j . Other attempts to generate TB-mixing relying on the non-abelian family 
symmetry A (27) give results which are also compatible with quark mixing and can 
in principle be embedded in a unified gauge theory [21J. However, the usefulness of 
parametrizations dictated by such symmetries is limited within the prescribed scenario 
and might not capture cases exhibiting other possible interesting properties beyond 
TB-mixing. Moreover, if we seriously wish to exploit the idea that some other underly- 
ing symmetry is found hidden behind the regularity of the neutrino data, pertubative 
investigations around the TB solution are highly unlikely to have a chance. We should 
not also ignore the variety of unified or string models where the existing vacua along 
flat directions usually break symmetries in a hard way and lead to complicated mass 
textures. Renormalization group effects, as well as instanton contributions [231 - 126] may 
further obscure the original symmetry. 

■^We note that the method developed here is general and can be applied equally well to the quark 
sector. 
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Taking into account the above considerations, we infer that the compUcated picture 
of the model building landscape as well as the sensitivity of the neutrino data on TB 
departures, call for a detailed and exact treatment of the neutrino sector. To accomplish 
this, in this section, we will develop a new formalism for describing 3x3 mass matrices 
and their corresponding diagonalizing transformations. In doing this, we will generalize 
a formalism which appears to be a special property of the original TB construction, 
namely the independence of the mixing angles from the eigenvalues. This will be 
a built-in property of our suggested formalism and will facilitate the analysis of the 
complicated structure of the leptonic sector. 

We wish to analyze general models based on GUTs, SUSY-GUTs and strings which 
predict a variety of fermion mass textures mf not necessarily symmetric or Hermitian. 
In the present analysis we will consider the Hermitian squares nifm}^ of the 3x3 fermion 
mass matrices which capture the physical properties of a whole class of fermion mass 
textures m/. 

A general Hermitean 3x3 matrix contains 9 independent elements and can be 
written as 

H = ihiU (2) 



where U a unitary matrix. Using the Caley-Hamilton theorem we can write 

H = bil + b2U + bsU'^ 



(3) 



where ^1,^25^3 are complex in general. Reversing the argument, we propose to write a 
Hermitean mass matrix M in the form 



M = bil + b2U + bsU^ 



(4) 



where U is a unitary matrix and, without loss of generality, we assume that det U = 1. 
Now, the standard CKM form for a unitary matrix contains four independent elements 
and has determinant one. Adding six degrees of freedom (d.o.f.) from the complex 6j 
coefficients we have a total of ten, so one d.o.f. is redundant, and can be removed by 
requiring one eigenvalue of U to be one as described in the next paragraph. 

Since U and M obviously commute, the above expression can be diagonalized by 
means of a similarity transformation. Thus, once we have expressed a given M in 
terms of J7, we can find its diagonalizing matrix simply by diagonalizing U. A diagonal 
unitary matrix is uniquely defined by 



lai 






















giaa 



(5) 



One phase can be absorbed into a redefinition of the coefficients 62 and 63 and taking 
into account the determinant condition, we end up with 














1 














(6) 
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where the ordering of the diagonal elements may vary. Thus, our unitary matrix can 
always be chosen to have one eigenvalue equal to one. Denoting by mi, 771-2, ?n.3 the 
(real) eigenvalues of M, we have the equations 



mi = 61 + b2e"^ + 636^'" 
m2 = bi+b2+ 63 
7713 = 61 + 626"*" + 636" 

The solution of the above system for 61 , &2 , &3 gives 



(7) 



-2ia 



bi 
b3 



1 2 

--CSC — 

4 2 



e 2'" 77li +62 7773 



2cosf 



7772 



1 
4 



,2 " 

2 



+- esc' — (e 



'(th-i - 1112) - (7772 - 7773)) 



a e 



--CSC — ■ 

4 2 



" 2'^(777l — ^-2) 



62" 



(7f72 - 7773) 



2cosf 



(8) 
(9) 
(10) 



Therefore, using this parametrization, we have succeeded to disentangle the mass eigen- 
values of M from the diagonalizing matrix. The eigenmasses 777j are given as functions 
of the coefficients bi and the phase a only. Consequently, for a given mass spectrum we 
may reconstruct the fermion mass texture by simply computing the coefficients bi from 
relations (|8l ll0p and a suitably chosen unitary matrix U. If for example, the mixing 
effects are accurately described by the experimental data, the mixing angles can be 
specified and U can be readily determined from the mixing matrix and the phase a. 
Next, we concentrate on the unitary matrix U assuming the standard parametrization 
in terms of three angles ^121^231^13 and a phase 5. We have 



U 



C12C13 

-C23S12 - Ci2Sl3S23e'^ 
-Cl2C23'Sl3e*'' + S12S23 



C13S12 
s 



C12C23 

-C23Sl2Sl3e 



e Sl2Sl3'523 
,5 



C12S23 



e-^'^si3 

Cl3'S23 
C13C23 



fill 



The requirement to have one eigenvalue equal to one leads to the constraint 

sin 6 sin ^12 sin ^13 sin ^23 = 0- 



(12) 



This condition is satisfied if one of the parameters vanishes generating four distinct 
structures for U : 



C12C13 

-C23'Sl2 - Cl2Sl3'S23 
-Cl2C23'Sl3 + S12S23 



Cl3 
-Sl3'S23e" 
-C23Sl3e*' 



C13S12 Sl3 

C12C23 - S12S13S23 C13S23 

-C23S12S13 - C12S23 C13C23 

e-^^si3 

C23 C13S23 1 ^12 

— S23 C13C23 



(13) 



(14) 
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Cl2 Si2 

-C23S12 C12C23 S23 
S12S23 -C12S23 C23 



^13 



C12C13 
-S12 
-Cl2Sl3e*' 



C13S12 
C12 

-Sl2Sl3e*' 



-i<5 



Sl3 





Cl3 



723 



(15) 



(16) 



In the present analysis, out of the four possible forms for U we choose the most ap- 
propriate cases and work out its implications on the leptonic sector. Since the mixing 
effects depend on the combined charged lepton and neutrino diagonalizing matrices, we 
treat them separately. 



2.1 The Neutrinos 



We start with the neutrino sector. Here we do not assume a specific embedding of the 
matrices in a given (GUT) model, thus the absolute scale of the neutrino eigenmasses 
is arbitrary and can be easily chosen in consistency with the /3/3-decay constraints. 

We will analyze the case where the unitary matrix associated to the neutrino mass 
texture Mi, is expressed as in ([4]) with U given by the specific form (fT4|l . i.e. 



U 



Cl3 
-Sl3'S23e* 
-C23Sl3e*' 





C23 
-S23 



e '■"Sis 
C13S23 
C13C23 



(17) 



By construction, the above matrix admits one eigenvalue equal to one, thus the eigen- 
values of U are 1 and e^*". The diagonalizing matrix for U is difficult to find in simple 
form, so we introduce the following parametrization. 



tan 9 



13 



2zi 



(18) 



tan 6*23 
5 



2ziZ2J{l + zl) {1-zl) 



+ 



vr 



z( + 2zfz^ 



(19) 
(20) 



Then, the diagonalizing matrix for U is 



1 

71 



pin 

Q 

P 



Z2 



iq^ Zl) 



{Z\Z2 
P 



-V2Z2 

V^qzi 



'"^[z2+iq' zi] 
-p {ziZ2 + ^) 



(211 
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where p, q are functions of 2:1^2 given by 



'1- J 



We can easily check that 

yJC/K = diagonal [e'", 1, e"'"] • (24) 
The eigenvalue a depends only on zi , Z2 

ia Zi - iZ2 _i Zi 

e = , or a = tan — 

Zl + iZ2 Z2 

The following relations are also useful: 

6*13 Zl cos ^ 

Zl = tan — , Z2 



2 ^Jzl+sin^ %i 



p=yi-zitan2^, g=ficot^. 

Using these equations, all the elements of ()2ip can be expressed in terms of the trigono- 
metric entries of the unitary matrix U. 



2.2 The Charged Leptons 



We now derive similar formulae for the charged leptons. As in the case of neutrinos, 
we choose to write the mass matrix in terms of a unitary matrix U in accordance to 
formula (jH). For reasons that will become clear later, we choose the ordering of the U 
matrix eigenvalues to be as follows: 



(25) 



If the eigenvalues of M are mi , m2 , 7713 the coefficients 6j in (jH) are given by equations 
analogous to (|8|9ll0p . 

For the case of charged leptons, we confine ourselves to orthogonal matrices. There- 
fore, out of the four possible forms for U we choose (jl3p 





"1 








u = 





















u 



C12C13 
-C23S12 - C12S13S23 
-C12C23S13 + S12S23 



C13S12 Si3 
C12C23 - S12S13S23 C13S23 
-C23S12S13 - C12S23 C13C23 



(26) 
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which corresponds to the structure Ui. Next, we use the fact that an orthogonal matrix 
can be written as H 

C/ = e""-^ = l + sinan-^ + (l-cosa) {n ■ 

where h = (ni,n2,n3) is a unit vector and the 3x3 matrices Si satisfy the conditions 



and are exphcitly given by 



si 



S2 



S3 



0' 

1 

-1 

■ 1 



-1 

1 0' 

-10 





The matrix U is diagonalized by means of the matrix 

1 



Vi (ni,n2,n3) 



V2Vnf + n| 
and we may check that 



2y + n^n^ n2n^ + ini 712723 — ini 



V^UVi = diagonal [l,e^",e" 



(27) 



(28) 



(29) 



(30) 



(31) 



2.3 The Leptonic Mixing Matrix 



In the previous sections we managed to obtain the diagonalizing matrices for the charged 
lepton and neutrino mass textures Vi and Vy in closed form. In accordance to standard 
notation, the leptonic mixing matrix is defined to be 



Vm = e'^V^V, 



(32) 



The phase factor is introduced in order to have the matrix determinant equal to one. A 
closer inspection reveals that all the elements of so derived Vm are complex. It can be 
shown [2^ however, that ()32p can be rendered equivalent to the standard form given in 



For a detailed mathematical analysis of the subsequent formalism in the context of the fermion 
meiss matrices see [14] . 
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(jlip . The proof goes as follows. A hermitean matrix M can by diagonalized by means 
of a unitary transformation 



= Diag [mi , m2 , ms] = D ■ Diag [mi , m2 , ms] • 
where -D is a miitary matrix of the form 



D = Diae 



(33) 



(34) 



This way, if [/ is a diagonalizing matrix so is UD . The lepton mixing matrix is given 

by 

Vm = V^V, (35) 
and taking the above into account, Vm can be equivalently written as 

Vm = D^V^V^C (36) 

where 

C = Umg[e'^\e'^\e'' 
D = Diag 



gidi gid2 ^id3 



If we require (p6]) be reduced to the standard form, the six phases can be uniquely 
determined. 

3 Analysis 

Using the above results, we will now proceed to determine possible deviations from the 
TB-mixing which fit the experimental data and determine the allowed range for the 
parameters zi,Z2,9,h in the neutrino and charged lepton sectors respectively. Before 
further pursuing the general case, we will first present the simplest and possibly the 
most elegant way of extending the Tri-Bi maximal mixing. 



3.1 Example: The minimal case 

We start with the neutrino sector and introduce values for the parameters zi^2 which 
are in accordance with the TB scenario. We put Z2 = —1 and we get tan ^23 = whilst 
for the eigenvalue of the unitary matrix U we get 



so that tan a 



and thus, 



a 



I — Zl 

-^13. This way the U matrix becomes 



U 



cos a 
1 

-ze~*^sinQ 



■ iO ■ 

-le sma 


cos a 



(37) 
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The neutrino mass matrix takes the simple form 



(mi + 7713) 


*^ (7713 — nil) 



ije 



(ma - mi) 





m2 
i (mi + ma) 



(38) 



For ^ = vr it is exactly the texture of the neutrino mass matrix introduced in the case 
of TB mixing [6] . The matrix (|38p is written in terms of its real eigenmass values with 
the {13}, {31} entries multiplied by a phase factor. 

Next, we proceed with the charged lepton sector. The TB-matrix [6] corresponds 

to 

We know however, that the TB case does not reproduce the experimental data since 
it predicts a zero ^13 angle. A minimal extension arises if we introduce in the neutrino 
mixing matrix the parameters 



Z2 



-1, 



7T + (f 



while keeping the charged lepton diagonalizing matrix as above. Then the mixing 
matrix is given by 



Vm 



ITT t<£_ -j- 

e 6 e 3 V, 



1 11 



' \V3' ^/3'^/3 

After some algebra and the removal of the redundant phase factors 
be brought into canonical form given by 



Vm 



2 ifi 

3COsf 



1 



2 ■ ip 

3 f 



sm 



sm 



\2^ q) V3 
2 &) 



I COS 



I COS 



(f+f) 

2 6/ 



(40) 

the matrix can 
(41) 



The experimental bounds are: 



0.0871557 < |sin 6*131 < 0.224931 
0.68728 < |tan6li2| < 0.713293 
0.213895 < |tan 6*231 < 1.09131 



and we have 



{Vm 


)ii 


{Vm 


)23 


{Vm 


)33 


{Vm 


)ll 


{Vm 


)l2 



sm 6^13 = 
tan 6*23 
tan O12 



2 Lf 

sm — 

3 2 



cos 



:f+f) 



cos (f - f ) 

1 



\/2 cos ¥ 



(42) 
(43) 
(44) 



(45) 
(46) 
(47) 
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Combining the above, we find that all the constraints are satisfied for 




(48) 



This is a rather interesting result: it states that TB-mixing can reconcile the neutrino 
data by a suitable choice of the phase parameter 9 parametrizing the neutrino diag- 
onalizing matrix. In the minimal case we are here dealing with, this phase coincides 
with the phase in the {13}, {31} elements of the neutrino mass matrix. In the original 
TB model this phase is simply taken to be = vr. When shifted by a value ip lying in 
the range (H8l) . neutrino data are exactly predicted. 

The advantage of this solution as compared to any perturbative approach around 
the TB-solution aiming to fit a nonzero 613 angle is rather obvious: indeed, it is shown 
that a nonzero ^13 angle that preserves the symmetric and zero form texture of the Mi 
and Ml, matrices can be naturally incorporated into the minimal TB-scheme. This is 
of crucial importance if we really wish to attribute their simple structure to some kind 
of discrete or other symmetry of the theory [9]- |21] . 

3.2 The general case 

We explore now regions of the parameter space which signal departures from the TB- 
case. Deviations can be easily obtained by assuming for example that 



in the charged leptons sector. Similarly, we choose to write ^2 = — l + fl'^ in the neutrino 
diagonalizing matrix (j2ip . In figure we plot the ranges of these parameters subject to 
the well known constraints of the mixing angles, while we keep 9 = tt. 

We observe that the allowed values of g lie in the range 0.05 < 5 < 0.6 , those of 
£ in the range —0.2 ^ e ^ 1 while acceptable values for zi cover a wider range lying 
from zi ~ 2.4 to large negative values zi ~ —10. It appears that there are wide regions 
in the parameter space consistent with data which significantly deviate from the TB 
mixing picture. 

Next, in order to determine mass textures related to possible exact symmetries, we 
scan the g, e ranges for fixed zi values. Here, we will concentrate in the subregions 
g ~ [0.2 — 0.6], e ~ [0.1 — 0.6] and search for values corresponding to exactly known 
trigonometric quantities. 

Let us choose e = l/^/2>■ This eliminates one entry in Vi which assumes the form 



Upon inspection, we observe that the remaining two parameters of the neutrino diag- 
onalizing matrix can be taken to be zi = tan ?f , 2^2 = — 75 • 




(49) 
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Figure 1: In this plot, deviations from TB-mixing are parametrized in terms of zi, 
= 1 + Z2 oi the neutrino and e of the corresponding charged lepton diagonalizing 
matrices, (see (I21|49j) .) 



With this choice, we can readily check that the matrix formed by the moduli of the 
elements of the leptonic mixing matrix is given by 



0.806056 
0.420639 
0.416337 



0.586939 
0.655601 
0.475067 



0.0759986 
0.627096 
0.775225 



pretty much close to the experimental data shown below 



0.546431 - 0.580416 0.03141 - 0.14091 ' 
0.63505 - 0.736914 



where the missing elements are determined by unitarity. For the above choice of param- 
eters, the charged lepton mass matrix obtained by substituting (|26p into the general 
form ([4]), is found to be 



Ml 



mi+rrt2+m3 


2mi— m2— m3 


„-m2-m3 


3 


3^2 


^ V6 


2m\—m2—ra^ 


im\+m2+mz 


•m2-m3 


3^2 


6 


' 2^3 


■m2-m3 


•m2-m3 






* 2^3 


2 



\ 



(50) 
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The matrix (jSOp exhibits an interesting structure and one could think of several ways 
to link it to a possible existence of underlying symmetries. For example, in cases of 
string derived models with several singlet fields (f),(p',... acquiring vevs, one defines 
expansion parameters e = {(j)) /M, e' = (</>') /M, . . . with M being the cutoff scale of the 
higher theory. Then, to leading order we get the hierarchy ^ m2 ^ mi and we 
can approximate this matrix by 

|e|2 ee' e 

Ml « ( ee' -|ef e' | m^ (51) 
e e' 1 

with e = i^^ ^,e' = and m^ a mass parameter related to charged lepton mass scale. 
The corresponding neutrino mass matrix is given by 



M^ 



I a+(mi+g_m2+m3) /j^ (mi - m3 ) +r7_ (mi - 2m2 +013 ) 4V2(mi-m3)-i(mi-2m2+m3) 

4 4 16 

mi+2m2+m3 »/3-(mi-m3)+7+(mi-2m2+m3) 

4 4 

a-(mi+g-).m2+m3) 
4 



where the dots stand for the corresponding complex conjugate entries and the various 
coefficients are 

6±V3 1±V3 1 I ^/2±^fl 32 ^ 

It is to be noted that all the off-diagonal entries of the neutrino mass matrix are 
expressed only in terms of the squared neutrino mass differences (note that we have 
assumed Hermitian squared mass matrices thus we have the correspondence mi -f-)- m^. ) . 
The resulting structure is now more complicated than the corresponding charged lepton 
one. This is of course to be anticipated in models employing the see-saw mechanism, 
since the effective neutrino mass matrix is a product of the Dirac and the heavy right 
handed Majorana neutrino mass matrices M^ oc m£)M^^m^. Depending on the specific 
structure of the hypothetical original theory, there are even more options to attribute 
this matrix to symmetry properties [28] . the analysis of this issue however goes beyond 
the scope of this letter. 



4 Conclusions 

In this letter we have investigated possible forms for the charged lepton and neutrino 
mass textures which can reconcile the experimental data on neutrino oscillations. In 
our analysis we have considered the hermitian squares of either mass matrix and used 
standard techniques to express each one of them as a second degree polynomial of 
a suitably chosen unitary matrix. Since the eigenmass dependence is encapsulated 
in the coefficients of this expansion only, we can express the neutrino mixing angles 
analytically, as functions of the parameters which define the unitary matrices that 
generate the charged lepton and neutrino mass textures respectively. Next, we may 
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use the available neutrino data on the mixing angles to constrain this parameter space. 
In particular, taking into account that the mass matrices suggested by Harrison ct 
al. are in good agreement with the Tri-Bi maximal neutrino mixing, we explored the 
parameter space for allowed deviations. We have found that the actual data including 
a non-vanishing ^13 angle can be nicely captured, by only introducing a single phase 
in the {13} and {31} entries of the neutrino mass texture in the original TB-scheme. 
Furthermore, upon varying the free parameters of our model in a wider range, we have 
found that neutrino data can be accommodated even for large deviations from the 
TB-matrices too. 
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